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ABST1ACT 

This  docnMat  describes  two  procedures  for  designing  recursive 
digital  filters  from  continuous- ties  filters  when  the  ratio  of  the 
e sapling  frequency  to  the  pole  frequency  is  saell.  The  coefficients  of 
the  proposed  digitel  filters,  which  are  derived  froa  the  step  and  rap 
invariance  of  the  corresponding  analog  filters,  have  been  determined 
for  real  and  complex  poles.  For  higher-order  filters  realised  in  e 
parallel  fora,  it  is  deaonstrated  that  the  dlscrete-tlae  transfer  func¬ 
tion  of  digitel  filters  obtained  by  the  step  and  raap  invariance  can  be 
derived  directly  froa  the  standard  s-transfomatlon  or  froa  the  pertlel 
fraction  expansion  of  the  continuous- tine  transfer  function.  The 
dlscrete-tlae  transfer  functions  of  the  step-  and  ranp-invsrlant  fil¬ 
ters  realised  in  e  cescade  forn  have  also  been  derived.  Finally,  the 
perforaance  of  these  aathods  is  daaonstrated  by  plotting  the  magnitude 
end  phase  responses  of  the  first-  and  second-order  digitel  filters. 

For  high-order  Butterworth  and  elliptic  filters,  the  magnitude  re¬ 
sponses  of  step-invariant  and  ramp-invariant  filters  ere  compered  with 
those  obtained  by  usual  aathods  such  as  the  standard  s  and  the  bilinear 
transformations. 

\  itOTMi 

Ce  document  dicrit  dawn  nftthodes  pear  concevoir  dee  f litres 
nwiriques  ricuralfs  h  partlr  des  f litres  anelogiques  lorsque  le  rap¬ 
port  entre  le  taux  d'dchantlllonaage  at  la  frequence  du  pile  est  fai- 
ble.  Las  coefficients  des  flltres  nsairlquee  proposis,  qui  provien- 
nent  de  1* Invariance  <  l'ichelon  at  <  la  range  des  f litres  anelogiques 
correspondents,  ant  %ti  ditemlnis  poor  les  phles  riels  et  complexes. 
Pour  lea  flltres  d'ordres  plus  Clevis  riellsCa  dens  an  riseau  paral¬ 
lels,  11  eat  dimontri  qua  la  fonctlon  de  tranefert  an  s  des  flltres 
nunirlques  cerectirisia  per  1' Invariance  i  l'ichelon  et  I  le  rampe  pent 
itre  diduite  dlreeteaent  de  la  transferals  ea  s  ou  de  la  decomposition 
en  fractions  pertielles  de  la  fonctlon  de  transfert  des  flltres  analo- 
glques.  On  a  aussl  didult  la  fonctlon  de  transfert  en  s  des  flltres 
invariants  I  l'ichelon  et  t  la  reaps  pour  one  realisation  dens  ua  r§- 
seau  an  sirie.  Final  ament,  le  performance  de  ces  nithodes  est  dfaon- 
trie  en  tragant  les  riponses  en  amplitude  et  ea  phase  des  flltres  nwai- 
riques  du  premier  et  du  denllae  ordre.  Four  des  flltres  Butterworth 
et  elllptlques  d'ordres  Clevis,  on  coopers  la  riponse  en  amplitude  des 
flltres  invariants  It  l'ichelon  et  1  la  ramps  avec  cells  obtenue  par  des 
nithodee  courentes  telles  qua  le  transformlc  an  s  at  la  transfomie 
blllnieire. 
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1.0  DtnOOPCTIOH 


As  with  analog  filters,  the  approximation  step  in  the  design  of 
digital  filters  is  the  process  whereby  a  realisable  transfer  function 
that  satisfies  prescribed  conditions  is  obtained.  Several  aethods 
permit  the  derivation  of  recursive  digital  filters  from  the  continuous** 
time  transfer  function  of  their  analog  counterpart  when  the  ratio  of  the 
sampling  frequency  to  the  pole  frequency  (fg/fp)  Is  sufficiently 
large.  The  methods  constitute  textbook  material  and  Include  the 
Impulse-Invariant,  matched  s  and  bilinear  transformations.  On  the  other 
hand,  they  suffer  from  serious  drawbacks  such  as  aliasing  and  frequency 
warping,  which  lead  to  Inappropriate  approximations  when  fg/fp  Is 
small. 


The  two  approximation  methods  described  in  this  report,  the 
step-invariant  and  ramp-invariant  transformations,  fill  this  gap  by 
providing  valid  matches  even  when  fg/fp  is  small.  The  proposed 
methods  are  an  extension  of  the  impulse-invariant  transformation,  filth 
this  technique,  which  is  also  called  the  standard  c- transformation,  the 
response  of  the  derived  digital  filter  to  an  Impulse  Is  identical  to 
that  of  the  sampled- impulse  response  of  the  continuous-time  filter.  In 
the  step* Invariant  transformation,  the  response  of  analog  filters  Is 
found  by  assuming  that  the  Input  is  approximated  by  a  sequence  of  steps 
whose  duration  is  set  to  the  sampling  period  and  whose  amplitude 
corresponds  to  the  instantaneous  value  of  the  input  signal  at  the  time 
of  sampling.  In  the  ramp- invariant  transformation,  the  Input  signal 
consists  of  a  sequence  of  ramps  that  Join  the  sampled  values.  This 
procedure  leads  to  the  determination  of  the  discrete- time  transfer 
function  of  digital  filters,  whose  response  is  Identical  to  that  of  the 
reference  analog  filter  In  relation  to  the  input  signal  (step  or  ramp). 
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In  thin  report,  discrete-time  transfer  functions  were  determined 
for  first-  and  second-order  filters  with  real  and  complex  poles.  Map¬ 
pings  from  the  continuous- time  (Laplace  transform)  and  discrete-time 
(s-transform)  transfer  functions  were  derived  for  the  step-  and  raap- 
inva riant  transformations.  The  methods  were  generalised  for  filters  of 
any  orders  with  poles  of  any  multiplicity  by  considering  partial  frac¬ 
tion  expansions  In  terns  of  first-  and  second-order  transfer  functions. 
This  approach  leads  directly  to  a  parallel  implementation.  The  discrete- 
time  transfer  function  of  the  step-invariant  and  ramp-invariant  filters 
realised  in  a  cascade  form  has  also  been  derived. 

This  theoretical  approach  was  validated  by  plotting  the  magnitude 
and  phase  response  of  first-  and  second-order  digital  filters  derived  fay 
the  step-  and  ramp-invariant  transformations.  These  responses  are  com¬ 
pared  with  those  of  the  analog  filters  and  the  impulse-invariant  digital 
filters.  Finally,  for  high-order  Butterworth  and  elliptic  filters,  the 
magnitude  responses  of  the  step-  and  ramp-invariant  filtera  are  compared 
with  those  obtained  by  the  impulse-invariant  transformation,  the  bilin¬ 
ear  tranaformatlon  and  the  original  analog  filter. 

This  work  was  performed  at  DREV  between  September  1982  and 
February  1983  under  PCN  21J05,  Guidance  and  Control  Concepts. 
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2.0  BACKGROUND 


Filtering  is  a  process  by  which  the  frequency  spectria  of  a 
signal  can  be  modified,  reshaped,  or  manipulated  according  to  some 
desired  specification.  It  may  imply  amplifying  or  attenuating  a  range 
of  frequency  components,  rejecting  or  isolating  one  specific  frequency 
component,  etc. 

The  digital  filter  is  used  to  process  discrete-time  or  sampled 
signals.  It  can  be  implemented  by  means  of  software  or  dedicated 
hardware,  and  it  can  be  represented  by  a  network  comprising  a  collec¬ 
tion  of  interconnected  elements.  The  analysis  of  a  digital  filter 
determines  the  response  of  the  filter  network  to  a  given  excitation. 
Its  design,  on  the  other  hand,  consists  in  synthetizing  and  implemen¬ 
ting  a  filter  network  so  that  a  set  of  prescribed  excitations  results 
in  a  set  of  desired  responses. 

pie  design  of  a  digital  filter  comprises  at  least  two  steps: 
the  approximation  and  realization  ones.  The  approximation  step  is  the 
process  of  generating  a  transfer  function  satisfying  a  set  of  specifi¬ 
cations  that  may  concern  the  amplitude,  phase,  and  possibly  time- ^ 
domain  response  of  the  filter.  The  realization  step  is  the  process  of 
converting  the  transfer  function  into  a  filter  network  using  inter¬ 
connected  unit  delays,  adders  and  multipliers. 

In  the  approximation  step,  the  wealth  of  knowledge  acquired  in 
the  design  of  analog  filters  may  be  transposed  to  the  design  of  their 
digital  counterparts.  The  discrete-time  transfer  function  (expressed 
in  z-transform)  is  derived  from  the  continuous-time  transfer  function 
(expressed  in  Laplace  transform)  of  a  reference  analog  filter  of  known 
characteristics.  Several  methods  are  currently  available  that  perform 
such  approximations:  the  Invariant- impulse,  matched  z-  and  bilinear 
transformations. 
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Once  the  discrete-tiae  transfer  function  of  the  digital  filter 
has  been  obtained,  this  transfer  function  can  be  converted  into  a  fil¬ 
ter  network.  The  realization  of  this  filter  network  will  not  necessar 
ily  operate  exactly  as  prescribed  because  the  previous  operations  were 
carried  out  asstaing  that  the  coaponents  to  be  used  are  of  infinite 
precision.  These  effects  need  to  be  studied  and  corrected  as  neces¬ 
sary.  The  realization  step  of  the  digital  filter  is  fully  described 
in  aany  textbooks  such  as  Refs.  1  and  2.  The  transfer  function  can  be 
broken  down  Into  slapler  transfer  functions  and  be  realized  directly. 
However,  if  the  transfer  function  is  given  as  a  sun  of  partial  frac¬ 
tions  or  as  a  product  of  first-  and  second-order  factors,  it  can  be 
realized  either  In  parallel  or  cascade  form  without  any  further 
■odi f lea t ions. 

2 .1  Impulse-Invariant  Transformation 

The  impulse- invariant  transformation  (Refs.  1  to  5)  yields  a 
digital  filter  with  an  Impulse  response  equal  to  the  sampled  impulse 
response  of  the  contlnous  filter.  In  this  transformation,  the  input 
signal  x(t)  is  sampled  at  frequency  fg.  The  response  of  the  analog 
filter  is  found  by  assuming  that  it  is  excited  by  a  sequence  of  im¬ 
pulses  equally  spaced  at  intervals  T  (the  sampling  period).  The  am¬ 
plitude  of  impulses  is  equal  to  that  of  the  sampled  values  of  x(t). 


Given  that  the  analog  filter  has  only  simple  poles,  its  transfe 
function  is  written  as: 


Ha(.)  -  A  +  l  — —  + 
A  °  i-1  *  +  pi 


N1 

l  { 

i-1 


-i—  +  } 

»  ♦  di  s  +  d* 


The  details  of  the  derivation  of  this  equation  can  be  found  in 
Appendix  A. 
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Hm  Inpnl  m  ’irngUat  transfornatioa  I(«)  of  eq.  1  |lm 


«*>  -  A  +  * 
0  i-1 


V 


l  -  o 


Ml 

“Pi*  -1  1»1 


N)i  ♦  Aii* 


-i 


1  -  (2«  cos  PjT)*  +  o 


-2^  ;  {*) 
-1  wi 1  —2 


where  s  -  T  bolng  tba  stapling  period 

*01  “  T  *11 
-«  T 

*1!  “  *•  [*!!  co*  ***  +  Mji  ,ln  *1*1 

*U  -  Cj  +  Ct  -  2  Bo(C1) 

*21  *  <C1  "  Cl>3  "  2  **«i> 

dl  “  *1  “  #1 

d*  -  conjugate  of  d^ 

The  aopplng  relation  between  the  ■  plane  and  the  s  plane  can  be 
deduced  directly  froa  eq.  2.  It  Is  written  as 


1 

•  +  Pt 


1 


T 


-  s 


[3] 


nils  aathod  yields  valid  astebes  with  the  corresponding  analog 
•7*tsn  only  at  high  stapling  rates  and  it  Is  satisfactory  only  when  the 
analog  systan  Is  sufficiently  bend  United. 
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2.2  Hite  bad  rftgjfoautoa 


The  matched  r-t ran* formation  la  a  technique  baaed  on  napping  the 
polea  and  seroa  of  tba  continuous-time  filter.  It  la  perforued  on  the 
eaaeade  fore  of  the  transfer  function  of  analog  flltera.  Aa  shown  in 
Appendix  A,  thla  fora  la  written  aa 


H(a) 


M' 

A  (‘  *  V 

M* 

A  <• +  pi> 


[*] 


The  discrete-time  transfer  function  of  eq.  4  is  given  in  Sef.  S  as 


M(«) 


(*  +  DL 


(*  -  e“"iT) 
(a  -  e‘plT) 


[5] 


where  L  la  an  Integer  whose  value  la  equal  to  the  nunber  of  aeroa  of 
H(s)  at  a  ■  In  this  transformation,  the  poles  of  digital  flltera 
are  Identical  to  those  obtained  by  the  lnpuls e-lnvarlant  tranafo na¬ 
tion.  However,  the  aeroa  do  not  correspond.  In  general,  the  use  of 
the  lupula e-lnvarlant  or  bilinear  trana formation  la  preferable  to  that 
of  the  matched  a-transforaatlon  (Kef.  1). 


2.3  Bilinear  Transformation 


A  bilinear  transformation  la  obtained  by  substituting  a  in  the 
continuous-time  transfer  function  H(s)  by 


a 


2  l_-  a 
T 

1  1  +  a 


-1 


-1 


[6] 


This  transformation,  which  is  used  to  circravent  the  aliasing 
problem  of  the  standard  r- transformation,  results  In  a  nonlinear  warping 
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of  the  fnqonej  scale  between  the  continuous-ties  frequency  and  the 
discrete- tine  one,  according  to  the  relation 


where  la  the  contlnuous-tlae  frequency  variable  and  la  the 
discrete-ties  one. 


3.0  8TEP-INVAEIANT  TRANSFORMATION 


The  step-invariant  trana forest Ion  la  an  extension  of  the 
Inpulse-lnvarlant  technique.  The  responae  of  the  analog  ayatan  la 
found  by  assuelng  that  the  Input  x(t)  la  approximated  by  a  sequence  of 
steps.  The  aeplltude  of  steps  corresponds  to  the  instantaneous  value 
of  x(t)  at  the  ties  of  sampling  and  their  duration  is  equal  to  saapllng 
period  T.  In  this  case,  the  input  signal  Is  written  as 

*<t)  -  J0  X<nff)  {|»(t  -  *T)  -  p(t  -  (n  +  1)  T)}  [8] 

where  u(t  -  nT)  ■  0  if  t  <  iff  and 
■  1  if  t  >  nT. 


In  the  following  sections,  an  invariant- step  transfornation  will  be 
derived  for  first-order  (real  poles)  and  second-order  (complex  poles) 
teres  of  the  partial  fraction  expansion  given  by  eq.  1. 

3.1  First-Order  Tares  (Real  Poles) 

Consider  a  real-partial  fraction  of  the  expansion  given  by  eq.  I 
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We  can  derive  free  eq.  9  e  difference  equation  that  will  give 
y[(n  +  l)Tl,  the  value  of  y(t)  at  tie*  t  .  and  froa  x(nT)  and 
y(nT),  the  values  of  x(t)  and  y(t)  at  tine  tQ.  The  transfer  function 
suet  he  rewritten  by  setting  the  Initial  condition  of  the  output 
variable  at  tiae  tQ.  Then,  the  filter  Is  stimulated  by  a  pulse  whose 
amplitude  and  duration  are  set  to  the  Instantaneous  value  of  the  Input 
signal  at  tiae  tn  and  to  the  sampling  period  respectively.  Substi¬ 
tuting  the  Input  variable  by  the  Laplace  transformation  of  a  pulse 
function  of  amplitude  x(nT)  and  multiplying  this  result  with  the  trans¬ 
fer  function  of  the  filter  yields  a  new  relationship.  Finding  its 
Inverse  Laplace  transformation  and  substituting  t  by  (n+l)T  gives  the 

value  of  the  output  at  tiae  t  .  , 

n  +  1. 

If  we  take  Into  account  the  initial  condition  of  the  variable 
y(t)  at  tiae  tQ(-nT) ,  noted  by  y(nT)  *(t  -  nT)  where 

6(t  -  nT)  -  1  If  t  -  nT 
•  0  elsewhere 


Thus,  the  transfer  function  of  eq.  9  is  transformed  as 


Y(m)  ,  MLtl  +  xCsT) .  .?.~nTa 

s  +  p  •  +  P 


Since  the  Input  variable  at  tiae  tQ  Is  a  pulse  of  amplitude 
x(nT),  Its  Laplace  transformation  Is 


e-(n  +  l)TSj 


[11] 


Substituting  eq.  11  In  eq.  10  and  finding  the  Inverse  Laplace 
transformation  of  this  relation  yields  the  following  time  function: 


x'i  *  :  ' 


a® 
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,<t)  -  (MSU<1  .  ,-K*  -  "))  ♦  r(d)  .-»“  *  *>}  ,<«  -  «i) 

* 

-  isMSi  (1  .  .-»<*  -  <•  ♦  W)  „(t  -  („  +  i)i)  !“ 

P 

Finally,  lat  ua  determine  tha  value  of  y(t)  at  tine  tft  +  ^  by 
replacing  t  by  (n  +  1)T  in  eq.  12*  This  laada  to 


y((n  +  1)T]  -  (1  -  a”1*)  +  y(nX)e"pT  [l3] 


Equation  13  la  a  finite  difference  equation  whose  discrete- tine 
transfer  function  is 


e/,\  m  IU1 
W  X(s) 


-*>  -1 


-  a-**."1 


By  conparlng  aqa*  14  and  9,  it  is  assn  that  tha  discrete- tias  transfer 
function  can  be  obtained  frou  tha  continuous- tine  tranafer  function  by 
using  tha  napping  relation 


ti  -  .-'V1 


a  +  p 


, -pT  -1 
1  -  a  r  a 


3.2  Second-Order  Te 


Nov  consider  one  coupler- partial  fraction  of  tha  expansion  given 


by  eq.  1 


a  +  d  a  +  d* 
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1 

*» 

t 


1 

I 

$ 


If  wa  apply  the  substitution  of  oq.  15  to  each  of  the  individual  te 
of  oq.  16,  no  find 


la-.-").-1  jo  -  r4*1).-1 

I-.-1.-"  *  l-.-1.-4** 


,  -2  .  _  -1 
V  +  v 


1  -  <2o“°T  coo  $  T)o"1  +  a“2aT*"2 
vhoro  *A  -  o”2aT  -  ^  o“®T  coo  6  T  -  Kj  a"®1  oln  6  T 


Kj  -  ^  -  Hx  o“®T  coo  S  T  ♦  *2  o“*T  oln  6  I 


Kj  -  (C/d)  +  <C*/d*) 


K2  -  <C*/d*  -  C/d)J 


Tho  otap-lavarlont  tranafotnatlon  of  analog  flltora  doocrlbod  by 
oq.  1,  vhooo  N  roal  polos  and  N1  conplox  poloo  ora  siapla,  lo  detar- 
adLnad  by  oobotltntlng  oqo.  IS  and  17  In  oq.  1.  This  givaa 

A1  "*iT  -i 

-i  (1  -  •  1  )■ 


.  r 


i-o*« 


_  —2  .  _  -1 
*A1*  *  Si* 

-a  T 

(2a  1  coo  6.T)o”1  +  o 


“2®iT  -2 
i  a 


l/> 


m\ 


¥ 


■ 


«rC,  i 


f$I 


M 


i.i 
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4.0  RAMP- INVARIANT  TRANSFORMATION 


The  e-transfora  of  this  technique  can  be  derived  fros  the  Laplace 
transform  by  applying  the  same  operations  that  yield  the  e-transformation 
of  the  step  invariance,  the  only  difference  being  that  the  analog  filter 
at  time  tQ  is  excited  by  a  ramp  such  as: 

*(tn)  “  {»[(n  -  1)T]  +  ^ - !2ll}  {w[t  -  (n  -  1)T]  ^ 

-  y(t  -  nT)} 

If  we  apply  the  following  procedure  to  a  real-partial  fraction 
(eq.  9)  of  the  expansion  given  fay  eq.  1: 

1)  Transformation  of  eq.  9  to  eq.  10:  this  operation  allows 

us  to  consider  initial  conditions  of  the  output  variable  at 

time  t  . 
n 

2)  Substitution  of  X(s)  in  eq.  10  by 

1)T1  +  x(nT)  -  x[ (n  -  1)t1 |  (e~(n  -  l)Ts_  e-nTs} 

•  Ts2 

which  is  the  Laplace  transformation  of  the  ramp  function 

(eq.  19)  at  time  t  . 

n 

3)  Solve  eq.  10  to  derive  a  relation  as  a  function  of  tine  and 
replace  t  by  off  to  get  the  value  of  the  output  at  time  tQ. 

A  difference  equation  that  gives  y(nT)  from  x(nT),  x[(n-l)T]  and 
y  [(n-l)Tl  can  be  derived.  Its  z-transformatlon  is 


for  liapit  poles,  tho  rap^lonrlMt  ttaufomtloa  of  analog 
flltor*  characterised  by  V  ml  polos  sol  A  eosplw  poioo  gloss: 


-»,r 


,<*)  •  *0  +  Jl  il 

♦  F— 


*1<1 '  »,t 


-)  +  (J 


-V 


-  V 


i  - 


Sa  ♦  ♦  «n«~* 


-o.T  ,  -2a  _  I  _ 

11  1  -  <2o  1  cos  i  t)*"1  ♦  o  1  a“2 


In  this  chap tar,  no  will  first  establish  s  procedure  that 
derives  the  step-  end  rasp-invariant  transformations  directly  from  the 
inpulse-invarlant  transformation.  Then,  thle  procedure  mill  be  used  to 
determine  the  step  and  rasp  Invariances  of  analog  filters  containing 
multiple  poles. 


The  napping  relation  between  the  s  and  s 
step  Invariance  is  from  eq.  15 


planes  that  gives 


s  +  p 


.  »/.)  » -  .~|IV1 

1  -  .-<Vl 


This  relation  can  be  rewritten  as 


•  ♦  P 


♦  <1 Ulrr- 


r~=r  - : .  • (1  -  *  >  z(«rr*> 


idiere  Z  is  defined  as  the  operator  corresponding  to  the  invariant- 
impulse  response  of  analog  systems  or  standard  (-transformation.  Also, 
the  mapping  relation  of  the  ramp  Invariance  (eq.  22)  can  be  expressed  as 


mt* 
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bffi- 
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Hass,  it  la  riam  that  the  iup  nd  mp  lwariaaeaa  can  1m 
larlaal  fees  tha  fafalaa  lrarlMca.  Uhls  epatatioa  eoaalata  la 
MltifiTiac  eh*  tMafar  function  of  at  wlag  filtar  Ka)  hp  ths 
laplaea  ttaaafonatiM  of  tha  aalt-atap  faaetiM  (i/s)  sc  tha  ssit-rae 
function  (1/s2),  finding  ths  stsaiirt  rtraaafatiatioa  of  this  talatlc 


finally,  disidisg  ths  result  bp  ths 


»- transformation  of 


ths  Mlt-st^  fssetiss  shieh  is  1/(1  -  s  )  #r  ths  unit-rasp  fssetiss, 
shleh  is  T  s"l/(l  -  s~1)2. 


Ths  ineariaot-step 
1(a)  of  ss  ssalst  spates 
pot  is  ths  foes 


sse  t(s)  ssd  ths  inuarisst-ranp  tm 
tcsssfsr  fssetiss  is  V(s)  ess  thso 


»(*)  •  {x£kl}]  a  -  «"l>  [W] 

ssd  i(s)  -  [x{Si}i}]  ft  :  j-1)2  [27] 

s*  h 

Ibis  procedure  siloes  ths  dsterslnstios  of  ths  stop  sod  rsap 
invariances  of  eaelog  filtsrs  Shoos  trsasfsr  function  contains  sultiple 
polos.  In  Appendix  A,  this  transfer  function  is  expressed  as: 


jp” 
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m# 
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1-1  *l  <•  +  p1>  1-1  <•  +  *tr  <•  +  dt> 


[2S] 


ahan  and  Ifl.^  in  Che  double  aiMatloa  indicetes  the  Multiplicity 
of  the  real  and  complex  poles. 


Iha  atandard  s-transforn  of  thaaa  functions  la  given  in  laf.  6 


as: 


(•  +  P  t) 


4+1 


4  "tpiT  -4 

T  *  W8 

(1  -  a'PiVV  +  1 


[29] 


where  D^a^)  is  defined  by 
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(4  -  1)1 
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(*  -  2)1 


ri  A 

where  a4  »  e  and  *  !•  For  coeplex  poles,  and  a^  are 


^i1 


replaced  by  d^  and  b^  respectively  (b^  -  a  ) 


The  derivation  of  the  step  invariance  can  be  found  by 
substituting  eq.  28  in  eq.  26.  This  yields 
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Chtpt*M  3.0  to  S.O  derived  tho  stop-  ud  rasp-invariant  digital 
filters  Croat  the  partial  fraction  expansion  of  the  contieuous-tlae 
transfer  enactions.  This  design  procedure  leads  directly  to  a  realisa¬ 
tion  la  a  parallel  fern.  If  the  cascade  or  direct  fore  la  desired.  It 
can  he  obtained  by  substituting  the  napping  relations  of  eqs.  15  and  22 
Into  the  transfer  function  of  analog  filters  expressed  as  a  cascade 
fore  of  first-  rad  second-order  blocks.  The  results  are  presented  In 
Appendix  B. 


6.0  CCHPA1I80W  OF  THE  METHODS 

The  Invariant-  lnpulse,  natchad  x-  and  bilinear  transfornatlone 
or  either  of  the  two  aethods  developed  In  this  report  lead  to  different 
discrete-tine  transfer  functions  and,  in  sene  eases,  to  different  real¬ 
isations.  At  sufficiently  high  fa/fp,  their  tine,  gain  and  phase 
responses  are  nearly  the  sane,  and  the  step-  and  ranp-invarlaat  trans¬ 
formations  will  be  unnecessarily  sophisticated.  In  this  chapter,  it 
•111  be  -down  that,  at  lew  values  of  f„/fp,  the  sophistication  inher¬ 
ent  in  these  neth-  ods  Is  necessary  for  transposing  the  character¬ 
istics  of  analog  filters  la  the  discrete-tins  donain. 

First,  the  approxiaation  of  the  first-  and  second-order  digital 
flltars  with  the  iapulse-,  step-  and  rasp-invariant  transfornatlone 
have  been  considered.  Tha  frequency  responses  In  eagnitude  and  in 
phase  of  the  transfer  functions  of  the  filters  obtained  by  eech  nsthod 
have  been  produced  for  the  second-order  band-pass  filter,  the  first- 
end  second-order  loir-pass  filters  and  the  flret-  and  second-order  hlgh- 
pess  filters.  The  frequency  responses  of  the  analog  filters  are  given 
as  a  reference  irtille  these  ef  the  digital  filters  are  given  for  a  set 
of  three  values  of  f§/fp:  100,  10  and  4.  In  all  cases,  the  pole 
frequency  fp  and  the  gain  of  the  analog  filter  at  fp  ears  set  to 
1  KBs  and  0  dB  respectively. 
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The  exercise  wee  repeated  for  higher-order  filters  such  as  the 
Butterworth  and  elliptic  low*pass,  high-pass,  band-pass  and  band-reject 
filters.  In  this  case,  only  the  frequency  responses  in  nagnitude  were 
produced.  However,  the  lap ulse-invar lent,  bilinear,  step-invariant  and 
rsmp-invarlant  transformations  were  coopered. 


6.1  Comparison  of  the  Methods  for  Deriving  First-  and  Second-Order 
Filters 

This  comparison  shows  the  differences  between  these  types  of 
digital  filters  in  the  time,  magnitude  and  phase  responses  when  small 
f,/fp  ratios  are  used.  It  leads  to  the  determination  of  the  best 
procedure  for  designing  high-order  filters  since  they  are  merely  a 
combination  of  first-  and  second-order  terms. 


First,  let  us  consider  a  second-order  band-pass  filter  with  the 
following  specifications:  quality  factor  Q  -  10,  resonant  frequency 
fQ  *  l  kHz  and  gain  at  resonant  frequency  -  1.  This  is  the  apecial 
case  of  eq.  A.5  where 


(«0/Q)a 


s2  +  (»0/Q)s  +  »Q2 


Figure  1  shows  the  performances  of  the  three  types  of  digital 

filters  for  various  sampling  frequencies  (f  ■  20  f  ,  10  f  ,  4  f  and 

s  o  o  o 

2.5  fQ)  when  the  input  is  excited  by  a  step  function  of  amplitude  equal 
to  1.  Each  time  the  ratio  f#/f  is  changed,  the  coefficients  of  the 
filters  are  recalculated  to  maintain  the  resonant  frequency  constant. 
From  Fig.  1,  ws  deduce  that  the  step-  and  ramp- invariant  procedures 
simulate  exactly  the  time  response  of  the  continuous  filter.  The 
impulse- Invariant  method  generates  an  offset  at  the  output  of  the 
filter  that  Increases  when  the  ssmpling  frequency  decreases. 


vv 
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Impulse-invariant  filter  Step-  (or  ramp-)  invariant  filter 


f  ■  1  kHz.  f  -  20  kHz 


f  *  1  kHz 


10  kHz 


f  -  4  kHz 


f  -  1  kHz,  f  -  2.5  kHz 

o  a 


FIGURE  1  -  Comparison  between  the  step-  (or  ramp-)  invariant  method 

and  the  impulse-invariant  method  in  the  time  response  to  a 
step-function  input  of  a  second-order  band-pass  filter 
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The  behaviour  of  the  tlae  riipooft  is  explained  by  examining  the 
magnitude  and  phase  responses  of  the  three  digital  filters  plotted  in 
Fig.  2  for  s  set  of  three  sampling  frequencies  (ffl  -  100,  10  and  4  fQ). 
The  responses  of  the  corresponding  analog  filters  are  also  plotted  in 
this  figure.  This  provides  a  reference  for  determining  the  performance 
of  the  digital  filters.  The  ramp- Invariant  method  produces  a  digital 
filter  with  a  frequency  response  in  amplitude  and  In  phase  that  corre¬ 
sponds  to  the  response  of  the  continuous  filter  for  a  frequency  band  of 
0  to  f#/2.  In  this  region,  its  frequency  response  does  not  depend  on 
the  sampling  frequency.  The  step-invariant  filter  has  a  magnitude 
response  similar  to  that  of  the  ramp-invariant  filter.  However,  its 
phase  response  does  not  match  the  desired  response  around  the  resonant 
frequency.  The  difference  between  the  analog  and  digital  responses 
does  not  exceed  30"  as  the  sampling  frequency  is  set  to  4  fQ.  The 
response  of  the  Impulse- Invariant  filter  is  identical  to  that  of  the 
analog  filter  only  in  the  frequency  band  located  around  the  resonant 
frequency.  Depending  on  the  sampling  frequency,  the  attenuation  of  the 
filter  in  the  frequency  hand  smaller  than  fQ  becomes  limited  to  a 
fixed  value  rather  than  following  the  slope  of  40  dB/decade.  For 
example,  at  fo/10,  the  digital  filter  gives  an  error  in  amplitude 
that  varies  from  12  to  18  dB  as  the  sampling  frequency  falls  between 
10  f0  and  4  fQ.  A  difference  between  the  digital  filter  and  its 
analog  counterpart  is  also  observed  in  the  phase  response.  It  can 
attain  90*  within  the  full  bandwidth  (f  /2)  when  f  is  set  to  4  f  . 


Legend  for  Figures  2  through  6 

.  Continuous- tins  filter 

Digital  filter  when  f#  ■  100  kHx 

—  — - Digital  filter  idien  f  •  10  kHz 

-  Digital  filter  when  f(  ■  4  kHz 
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FIGURE  3  -  Magnitude  and  phase  responses  of  a  first-order  low-pass 
digital  filter  obtained  by  the  impulse-,  step-  and  ramp- 
invariant  methods  for  various  sampling  frequencies 
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FIGURE  4  -  Magnitude  and  phase  responses  of  a  second-order  low-pass 
digital  filter  obtained  by  the  Impulse-,  step-  and  r amp- 
invariant  methods  for  various  sampling  frequencies 
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FIGURE  5  - 


Magnitude  and  phase  responses  of  a  first-order  high-pass 
digital  filter  obtained  by  the  impulse-,  step-  and  ramp- 
invariant  methods  for  various  sampling  frequencies 


Frequency  »•»  H*  Frequency  ifi  Mt 

Magnitude  and  phase  responses  of  a  second-order  high-pass 
digital  filter  obtained  by  the  inpulse- ,  step-  and  ramp- 
invariant  methods  for  various  sampling  frequencies 
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The  discrepancies  In  Che  nsgnlfcude  end  phase  responses  also 
exist  for  Che  low-  and  high- pass  filters.  Figures  3  and  4  Illustrate 
these  responses  of  the  first-  and  second-order  low-pass  filters  for  the 
step-,  reap -  and  Inpulse- Invariant  procedures  when  the  s sapling  frequen¬ 
cies  are  set  to  100,  10  end  4  fp.  The  pole  frequency  end  the  gain  of 
the  filters  are  neintalned  cons tent  at  1  kHz  and  1  respectively.  The 
quality  factor  Q  of  the  second-order  filters  Is  set  to  10.  The  Inpulse- 
invarlent  filter  generates  en  error  in  asgnltude  that  varies  with  the 

sanpllng  rate.  When  f  is  set  at  4  f  ,  the  error  Is  United  to  6  dB  for 

s  p 

the  first-order  filter  and  to  3  dB  for  the  second-order  filter  in  the 
frequency  hand  of  0  to  tJ2»  On  the  other  hand,  the  low-pess  digital 
filters  fron  the  step-invariant  and  ranp-lnvarlant  nethods  slnulate 
exactly  the  characteristics  of  their  analog  counterparts  for  a  frequency 
band  lower  than  fg/2,  except  that  the  phase  of  the  second-order  step- 
invariant  filter  is  slightly  disturbed  around  the  resonant  frequency. 

Finally,  the  nagnltude  and  phase  response  of  first-  and  second- 

order  high-pass  filters  are  illustrated  la  Figs.  5  and  6  for  the  sees 

operating  conditions  as  for  the  low- pass  ones.  For  the  first-order 

uodel,  the  preservation  of  the  analog  filter's  charecterlstlcs  by  the 

lupulse- invariant  net hod  requires  a  high  sanpllng  rate  (>  100  f  ).  For 

P 

the  second-order  uodel,  the  analog  filter's  characteristics  are  pre¬ 
served  on  a  narrow  frequency  band.  The  attenuation  of  the  amplitude 
becoaes  Halted  to  e  constant  value  in  the  frequency  band  lower  than 
f(/2  rather  than  following  a  constant  slope  of  40  dB/decade.  The  atten¬ 
uation  value  and  the  starting  frequency  of  this  phenomenon  vary  with  the 
sanpllng  rate.  The  step  invariance  preserves  the  characteristics  of  the 
first-order  high-pass  filter  within  the  frequency  band  of  0  to  tjl  but, 
at  snail  f ^/f ^  ratios,  a  constant  shift  in  the  nagnltude  response  is 
obaerved.  This  shift  varies  in  function  of  the  sanpllng  rate  and,  for 
example,  it  is  about  8  dB  at  f(  ■  4  fp.  In  addition,  an  error  in  the 
phase  response  exists  around  the  pole  frequency.  For  the  second-order 
filter,  the  constant  amplitude  attenuation  encountered  in  the  impulse- 
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Invariant  filter  la  replaced  b y  an  atteaaatloa  reduced  eo  e  slope  of 

20  dB/octava.  Hie  filter  also  generates  a  phase  error  of  90*  la  the 

useful  frequency  bend  for  stapling  rates  aaaller  than  10  f  .  The 

P 

aagnltude  and  phase  responses  of  the  first-  and  second-order  high-pass 
digital  filters  free  the  reap  Invariance  correspond  to  the  responses  of 
the  high- pass  analog  filter,  provided  that  the  frequency  band  Is  Ha¬ 
lted  to  f#/2. 

Thus,  step-  and  rasp- Invariant  filters  give  aore  accurate  fre¬ 
quency  response  in  aagnltude  and  In  phase  than  the  lapuls e-lnvarlant 
filter  for  all  the  cases  of  first-  and  second-order  filters  when  saall 

f  /f  ratios  are  used.  In  this  condition,  the  iapulse  invariance  can 
•  P 

be  perforaad  only  on  the  second-order  low- pees  filter.  The  step  In¬ 
variance  gives  accurate  frequency  responses  for  low- pass  and  band-pass 
filters  with  saall  phaae  errors  around  the  coaplax  poles.  Finally,  the 
reap  Invariance  produces  a  close  natch  between  the  frequency  reaponses 
of  digital  and  analog  filters. 

6.2  Conparlson  of  the  Methods  for  Deriving  Higher-Order  Filters 

The  next  set  of  figures  (Flga.  7  to  15)  shows  the  frequency 
aagnltude  characteristics  of  digital  filters  obtained  by  the  step-, 
reap-  and  iapulse- invariant  and  the  bilinear  transforaations  of  high- 
order  Butterworth  and  elliptic  filters.  For  further  details  about 
Butterworth  and  elliptic  filter  synthesis,  see  Ref.  7.  In  all  exaa- 
ples,  the  response  of  digital  filters  Is  presented  for  various  sanpling 
frequencies  (1000,  100,  10  and  4  kHz)  and  It  is  accoapanled  by  the 
response  of' the  corresponding  analog  filter.  Figures  7  to  10  Illus¬ 
trate  fifth-order  low-pass  and  high- pass  Butterworth  and  elliptic 
filters.  The  cut-off  frequency  and  the  gain  of  these  filters  are 
aalntalned  at  1  kHz  and  1  respectively.  In  addition,  for  the  elliptic 
filters,  the  stop-band  attenuation  is  at  least  40  dB  and  the  pass-band 
ripple  Is  llaltsd  to  3  dB.  The  band-pass  and  band-reject  Butterworth 
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sad  elliptic  filters  of  Pits.  11  to  14  were  derived  from  s  transforna- 
tioo  perforaed  on  the  corresponding  low pees  filters  (tef.  7).  the 
lever*  sad  upper-cutoff  frequencies  of  these  flltere  ere  located  at  700 
8s  sad  1*4  Ub.  the  polee  sad  seros  of  theae  analog  filters  ere  giden 
la  dppendlx  C.  the  cutoff  frequency  of  the  high-  sad  low- pass  filters 
am d  the  center  frequency  of  the  hand-pass  and  hand-reject  ones  ere 
noted  tjr  fp.  In  addition,  the  step  and  r amp  Invariances  of  analog 
flltere  preeented  in  this  group  of  figures  were  perforaed  on  the 
partial- fraction  expansion  which  yielded  a  parallel  realisation. 

this  set  of  figures  shows  that  the  lapulse-invarlant  transfor¬ 
mation  can  be  applied  only  to  Buttervorth  lowpass  and  band-pass 
filters.  In  the  flrat  case,  a  close  aetch  between  the  analog  and 
digital  filter  is  obtained.  In  the  aecond  case,  the  uagnltude  responee 
in  the  low- frequency  band  does  not  correspond  to  that  of  its  analog 

counterpart  when  the  ratio  f  It  becoaes  saaller  than  10.  In  all  the 

•  P 

other  esaaples,  the  Impulse-Invariant  aethod  is  unacceptable.  For 
elliptic  low-pass  and  hand-pass  filters,  the  equlripple  character  of 
the  stop-band  response  was  destroyed  by  the  aliasing  effect.  Finally, 
the  aliasing  renders  this  transformation  entirely  useless  as  a  digital 
band-reject  or  high-pass  filter  (Butterworth  or  elliptic). 

The  bllinearly  transformed  filters  are  essentially  Identical  to 
the  original  analog  filter  dien  the  ratio  fg/fp  is  greater  then  10. 
However,  the  polee  of  the  digital  filters  are  aoved  in  relation  to  the 
s sapling  frequency.  This  effect,  called  nonlinear  frequency  warping, 

is  usually  observed  when  the  s sapling  frequency  is  saaller  than  10  f  . 

P 

Thus,  this  aethod  does  not  correctly  preserve  the  magnitude  response  of 
analog  filters  at  snail  values  of  fs/fp. 

The  step  and  reap  invariance  of  lowpass  and  band-pass  filters 
produce  digital  filters  whose  aagnltude  response  does  not  depend  on  the 
sampling  frequency  la  the  frequency  band  of  0  to  f #/2.  For  high-pass 


and  hrthwjirt  f Hears,  tha  angaltude  mpem  of  the  step-invariant 
IHtttl  f  11  tars  la  unacceptable  at  mil  valuu  of  tia  ratio  f#/fp. 
ta  rafrlwartaat  method  can  bo  aaad  to  digitise  hl^t-pui  and  biad* 
reject  flltoro  but  tha  magnitude  roapoaaos  of  analog  and  digital 
system  am  not  idontlcal  ohm  tha  ratio  f|/fp  becomes  mailer  than 
10*  la  thla  case,  tha  equi ripple  eharoctor  of  elliptic  high- pass  and 
band- reject  flltoro  la  not  correctly  simulated.  Iho  aagaitado  roaponao 
of  tha  lutterworth  high-pass  filter  in  tha  loo- frequency  band  is  not 
preserved  by  the  ramp-invariant  filter.  Moreover,  the  attenuation  of 
tha  bead-reject  Bettenorth  filter  in  the  stop-hand  region  is  United 
to  a  value  that  varies  as  a  function  of  tha  saapllng  frequency. 

If  the  realisation  in  cascade  is  used  for  the  rap  Invariance  of 

high-pass  and  band-reject  filters,  bettor  approximation  of  the  analog 

filters  is  obtained.  Figure  IS  reveals  that  the  digital  lutterworth 

filters  preserve  the  aagnitude  response  of  their  analog  counterparts. 

However,  a  snail  shift  in  the  aagnitude  roaponao  of  the  digital  high- 

pass  filter  can  be  observed  when  the  ratio  f/f  becomes  smaller  than 

•  P 

10.  For  the  elliptic  high-pass  and  hand-reject  filters,  a  close  natch 

between  analog  mad  digital  eystaas  exists  if  the  ratio  f  /f  is  aain- 

•  P 

talnad  above  10.  Cascade  realisation  can  also  be  performed  on  low-pass 
and  band-pass  step-invariant  filters  but  the  aagnitude  response  of 
these  filters  is  less  accurate  than  that  obtained  by  the  rsap-lnvarlant 
filters  realised  in  a  parallel  fora. 


Legend  for  Figures  7  through  15 
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Digital  filter  when  f 
Digital  filter  when  f 
Digital  filter  when  f 
Digital  filter  when  f 
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-  10  kHs 

-  4  kHs 


liptlc  low-pass  digital  filter 
a riant,  the  bilinear,  the  step- 
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Oti* 


FIGURE  9  -  Magnitude  response  of  e  Buttemorth  high- pass  digital 
filter  obtained  by  the  lepulse- Invariant,  the  bilinear, 
the  step-  and  the  map-invariant  transforestions  for 
various  sampling  frequencies 


i  response  < 
by  the  lspi 
rsap-lmrarli 


FIGURE  II  -  Magnitude  response  of  s  Butterworth  bend-pass  digital 

filter  obtained  by  the  laqpulse-invsrlant ,  the  bilinear,  the 
step-  and  the  raap-invarlant  transformations  for  various 
sampling  frequencies 
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7.0  OOWCLPSIOM 


Two  aethods  for  approximating  digital  r« curs It*  filters  fro* 
analog  filter*  have  been  described.  The  invariant-step  response  is 
maintained  in  one  case  while  the  invariant-rasp  response  is  used  in  the 
other  case.  The  discrete-tine  transfer  fractions  of  digital  filters 
obtalnsd  by  these  aethods  were  derived  for  parallel  and  cascade 
realisations. 

The  perforaaace  of  the  proposed  aethods  was  first  detemlned  by 
giving  the  aagnitude  and  phase  responses  of  the  first-  and  second-order 
filters.  It  is  shown  that,  for  low-pass  and  band-pass  filters,  the 
step-invariant  procedure  gives  accurate  aagnitude  responses  that  do  not 
vary  with  the  saapllng  frequency  in  e  frequency  band  of  0  to  f^/2.  The 
reap  invariance  gives  aagnitude  responses  as  accurate  as  those  produced 
by  the  step-invariant  asthod  for  low-pass  and  band-pass  filters  and  it 

can  be  used  to  reproduce  perfectly  highrpass  filters  for  ratios  f  /f 

•  P 

saaller  than  10.  In  addition,  a  close  astch  in  the  phase  response 
between  the  rasp-invariant  digital  filter  and  its  analog  counterpart 
can  be  observed  for  all  cases  of  first-  and  second-order  filters. 

Finally,  a  comparison  between  the  step-,  ramp-  and  impulse- 
invariant  and  the  bilinear  transformations  was  performed  on  high-order 
Butterworth  and  elliptic  filters.  At  smell  sampling  rates,  the 
Impulse-invariant  procedure  realises  s  close  astch  with  the  analog 
filter  in  the  aagnitude  response  only  for  the  Butterworth  low-pass 
filter.  In  the  bilinear  transformation,  the  poles  of  digital  filters 
are  moved  in  relation  to  the  saapllng  frequency.  This  distortion, 
caused  by  nonlinear  frequency  warping,  bscoaas  significant  when  the 
ratio  tjt p  becomes  smaller  than  10. 
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The  step  sad  rasp  invariance  of  low-pass  and  band-pass  filtara 
in  a  parallel  implementation  producer  digital  filters  whose  magnitude 
response  does  not  depend  on  the  sampling  frequency  in  the  frequency 
of  0  to  fs/2*  The  magnitude  response  of  the  step-invariant  high- 
pass  and  hand-reject  filters  realised  in  a  parallel  or  a  cascade  form 
is  unacceptable  for  snail  ratios  f#/f  (<10).  The  ramp  Invariance 
realised  in  a  cascade  form  can  be  used  to  digitise  high-pass  and  hand- 


reject  filters.  The  ratio  f  /f  can  be  decreased  to  10  for  Butterworth 

•  P 

filters  but  the  equirlpple  character  in  the  stop-band  response  of 


elliptic  filters  is  preserved  by  a  ratio  fs/fp  of  at  least  10. 


The  step-  and  ramp-invariant  methods  are  thus  less  sensitive  to 

frequency  folding  in  comparison  with  the  impulse-invariant  method. 

Furthermore,  step  invariance  can  be  applied  to  transfer  functions  in 

which  the  degree  of  the  denominator  must  exceed  that  of  the  maeerator 

by  at  least  one.  In  the  rasp  invariance,  the  numerator  degree  can  be 

as  high  as  the  denominator  degree.  In  addition,  when  the  ramp- 

invariant  transformation  is  Implemented  in  a  parallel  form  for  the 

low-pass  sad  hand-pass  filters  and  in  a  cascade  form  for  the  high-pass 

and  band-reject  filters,  the  resulting  natch  becomes  better  than  the 

one  with  the  bilinear  transformation  for  ratios  f  /f  smeller  than  10. 

•  P 


/  / 
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APPENDIX  A 


Representation  of  the  Continuous-Tine 


Transfer  Function  of  Anal 


Let  us  consider  that 


i  -  Ill!  -  X  b*" 
VB)  x(.)  y  vi 


[A.l] 


la  the  transfer  function  of  a  realisable  analog  filter  If  coefficients 
a^  and  are  all  real,  and  the  degree  of  the  nunerator  Is  snaller 
than  or  equal  to  the  degree  of  the  denominator. 

Equation  A.l  can  be  expanded  into  partial  fractions  and,  after 
separating  the  real  from  the  cosplex  poles.  It  yields: 


U1  Alk  t  P1  ^1  B2ik*  *  Bllk 

1  +  p^k  1*1  *1  (,2  +  2ats  +  «2  +  S2)k 

y  f1  V.  +  f  21  — [A.3] 
i“l  *1  (S  +  Pl)k  (s  +  d^k  (s  +  d*)k 


H(s)  -  AQ  + 


H(s)  -  Aq  + 


Wher*  Clk  "  (B2ik/2^  +  J  ^B2ik“l  "  *11^ /2  Bl^ 


di  -  at  -  j  st 


The  *  Indicates  the  coaqplex  conjugate.  The  nultlpllcity  and  Ml^ 
of  real  poles  p^  and  complex  poles  (d^,  d*>  are  related  to  the  degree 
of  the  denominator  of  H(s)  by 
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*'  ■  £  mi + “t 

If  va  U8um  that  all  raal  and  coaplex  poles  are  slaple,  eq.  A.3 
becoaea  slaply 


"<•>  -  *o  +  i,  7 


N1 

+  y 


[A 


Equation  A. 5  covers  low-pass,  high-pass,  band-pass  and  band- 
reject  filters  aaong  the  well- known  filter  classes  such  as  the 
Butterworth,  Bessel,  Chebfshev  and  elliptic  filters. 


Equations  A.2  and  A.3  are  suitable  for  representing  parallel 
realisation  of  analog  filters.  For  serial  or  cascade  realisation, 
eq.  A.l  aust  be  expressed  as 


H(s) 


i5i 


N' 

1-1 


(•  + 

(•  ♦  Pt) 


[A 


If  ue  separate  the  first-  and  second-order  terns,  eq.  A.6  bee ones 

2 


.  .  ,  S*  *21*  *  *11,  ™  ’si*  +  *«•  +  *31 

H(s)  -  K  in  - - - 1  E  —s - s - z 

•  *-»  *  +  ',l  »*»  .2  +  2«,.  +  a  *  +  «2 


[A 


where  N*  -  NA  +  2NB. 


Equation  A. 7  can  also  be  written  as 
NA 


ci*i 


«•>  •  *.  [,iT  <*21 ♦  )]  ill  t»5i +  rra; +  7%)  1* 
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APPENDIX  B 


Realisation  of  the  Step-  and  Ranp-Invarlant 


_ „ital  Filters  In  a  Cascade  Forn 

If  the  cascade  or  direct  fora  la  desired.  It  can  be  obtained  by 
substituting  the  napping  relations  of  eqs.  15  and  22  Into  the  transfer 
functions  of  analog  filters  expressed  In  cascade  forn  (eq.  A. 8). 


For  such  a  realisation,  the  step-invariant  transfornatlon 


becomes 


.  HA 

S(x)  -  [£  Aj,  + 


_p  IJ 

{<AU/Pl)  -  AjJ  (1  -  e  1  s'1) 

i  ~P*T  -1 

1  -  e  s 


V  “2  V  -1 

All*  +  KB11* 


cos  BjTs”1  +  e 


.  J*  -  .  *C11  *  *D11*  *  cllZ 

51  -a.T  i  “2a. T  , 

1  -  2e  1  cos  B1Ts"1  +  e  1  z~£ 


[B.2] 


where  K^,  end  ere  directly  obtelned  by  substituting  C  by  In 
eq#  23# 


-  nwInVVX-V-.-’v’  M- V' 


.  AfvfeM’V  /vi  '*rin 


-.y>g  -» •*  **.*>#:*»*": 


fi,~  V-  k-”.Nr' 
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APPENDIX  C 


Continuous-Time  Transfer  Functions  of 


tbs  Buttercorth  and  Elliptic  Filters 


This  appendix  gives  the  continuous- time  transfer  function  of  the 
Buttexvorth  and  elliptic  filters  used  in  Chapter  8.0. 


The  transfer  function  of  these  filters  is  written  as 


H(s)  -  K 


lit  +  ■!> 

•  r  <•  +  pi> 

i«i 


[c.i] 


1)  Fifth-order  Buttervorth  low-pass  filter: 

-  cutoff  frequency  ”  1000  Hz 

-  M'  ■  rnnber  of  zeros  -  0 

-  N'  ■  nuaber  of  poles  ■  5 

-  Kg  -  9.79  x  1018 

-  p1>2  -  1941  ±  j  5975 
P3’A  -  5083  ±  j  3693 
p5’«  6283 

2)  Butterworth  high-pass  filter  derived  from  the  low-pass 
filter: 

-  cutoff  frequency  ■  1000  Hz 

-  M'  -  5,  N’  -  5,  K#  -  1 

-  poles  are  identical  to  those  of  the  low-pass  filter 

-  zeros  are  all  located  at  0 


m&m, 


*-"V,W,VA 

.scn'o:-. 


w- 

v!k 
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3)  Buttervorth  band-pass  filter  derived  from  the  low-peas 

filter: 

-  lower-cutoff  frequency  ■  700  He 

-  upper-cutoff  frequency  “  1400  Hz 

-  M*  -  5,  N*  -  10,  -  1.731  x  1018 

-  zeros  ere  all  located  to  0 

-  Pl  2  -  2221  i  J  5877 
P3’4  »  1416  ±  J  4866 
p5’6  -  2177  ±  j  7478 
p7’g  .  466.6  ±  j  4484 
P9’l0  "  906.3  ±  j  8709 

4)  Buttervorth  band-reject  filter  derived  fron  the  low-pass 

filter: 

-  lower-cutoff  frequency  •  700  Hz 

-  upper-cutoff  frequency  •  1400  Hz 

-  M'  -  10,  r  -  10,  Kg  •  1 

-  poles  are  Identical  to  those  of  the  band-pass  filter 

-  10  zeros  are  located  at  0  ±  J  6283 


5)  Fifth-order  elliptic  low-pass  filter: 


-  cutoff  frequency  ■  1000  Hz 

-  M'  -  4,  N'  -  5,  K8  -  240 

-  ■1  2  -  0  ±  j  9309,  «3  4  -  0  ±  J  6865 

-  p1’2  -  775.2  ±  J  4367,’  P3  4  -  155.5  ±  J  5802 
P.’-  1479 


Km 
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6)  Elliptic  high-pass  filter  derived  f roe  the  low-pass  filter: 

-  cutoff  frequency  -  1000  Hz 

-  M*  -  5,  N»  -  5,  K,  -  1 

"  *1  2  "  °  1  J  4240*  "3  4  “  0  *  -J  5750 

•s’"  ° 

-  Px  2  -  182.3  ±  J  6799,  P3  4  -  1555  ±  j  8763 
P5’-  26689 

7)  Elliptic  band-pass  filter  derived  f roe  the  low-pass  filter: 

-  lower- cutoff  frequency  ■  700  Hz 

-  upper-cutoff  frequency  “  1400  Hz 

-  M'  -  9,  N'  -  10,  -  169.72 

-  *1  2  "  0  1  J  10324,  *3  4  '  0  4  j  3801 

■j*4  -  0  ±  j  9163,  mj  ’  -  0  ±  J  4308 

■9’-  ° 

-  pl  2  “  379  ±  J  4557,  P3  4  -  72.08  t  j  8660 

P5’6  -  208.6  ±  j  4920,  p?  g  -  339.5  ±  J  8008 

P9*10  -  5229  ±  J  6261 

8)  Elliptic  band-reject  filter  derived  froe  the  low-pass 
filter: 


-  lower-cutoff  frequency  ■  700  Hz 

-  upper-cutoff  frequency  ■  1,400  Hz 

-  M’  -  10,  N*  -  10,  %  -  1 

-  *3  2  -  0  ±  J  6283,  «3  4  -  0  ±  j  7959 
■5*6  -  0  ±  j  4960,  By*  -  0  ±  J  8636 
<10  -  *  *  J  4570 

-  px  -  16,476,  p2  -  2396 

P3  4  -  41.4  ±  j  4323,  p5  6  -  87.4  ±  J  9131 
P7>g  -  306.1  ±  j  3889,  p9  1Q  -  793.7  ±  j  10086 


